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We continue our studies Q of an one-dimensional anisotropic exclusion process with parallel dy- 
namics describing particles moving to the right on a chain of L sites. Instead of considering periodic 
boundary conditions with a defect as in m we study open boundary conditions with injection of 
particles with rate a at the origin and absorption of particles with rate (3 at the boundary. We 
construct the steady state and compute the density profile as a function of a and (3. In the large 
L limit we find a high density phase (a > (3) and a low density phase (a < (3). In both phases the 
density distribution along the chain approaches its respective constant bulk value exponentially on 
a length scale £. They are separated by a phase transition line where £ diverges and where the 
density increases linearly with the distance from the origin. Furthermore we present exact expres- 
sions for all equal-time n-point density correlation functions and for the time-dependent two-point 
function in the steady state. We compare our results with predictions from local dynamical scaling 
and discuss some conjectures for other exclusion models. 
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I. INTRODUCTION 



We study an one- dimensional totally asymmetric exclusion model where particles are 
injected stochastically at the origin of a chain of L sites, move to the right according to 
rules defined below and are removed at its end, again according to stochastic rules. Each 
site of the chain can be occupied by at most one particle. Among the interesting features 
of such exclusion models is the occurence of various types of phase transitions which arise 
from the interplay of the bulk dynamics with the boundary conditions [|I|]- || and their close 
relationship to vertex models M, growth models J7J, and, in the continuum limit, to the 
KPZ equation || and the noisy Burger's equation. 

Exclusion models can be divided into four classes according to the dynamics (parallel or 
sequential) and the boundary conditions (periodic with conservation of the number of parti- 
cles and (possibly) a defect or open with injection and absorption of particles). According to 
this classification we call them p/p- models (parallel, periodic), p/o-models (parallel, open), 
s/p-models (sequential, periodic) and s/o-models (sequential, open). A s/p-model with a 
defect has been studied numerically by Janowsky and Lebowitz ||, without a defect it was 
solved by a Bethe ansatz by Gwa and Spohn ||. The s/o- model was studied numerically by 
Krug 0, later the exact solution with the full phase diagram was found In addition 

to the density profile all equal-time n-point density correlation functions in the steady state 
were determined IfLOfl . Previously we solved a p/p- model with a defect [||] with Bethe ansatz 
methods and obtained the density profile and the equal-time two-point correlation function 
in the steady state. 

Here we discuss the p/o- model with the same parallel dynamics as in [||] but with open 
boundary conditions where particles are injected at the origin with a rate a and are removed 
at the (right) boundary with rate (3. The bulk dynamics of our model are deterministic and 
defined as follows: Each site x on the ring (1 < x < L) is either occupied (r x (t) = 1) or 
empty (r x (t) = 0) at time t. The time evolution consists of two half time steps. In the first 
half step we divide the chain with L sites (L even) into pairs of sites (2,3), (4,5), . . ., (L, 1). 
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If both sites in a pair are occupied or empty or if site 2x is empty and site 2x + 1 occupied, 
they remain so at the intermediate time if — t + 1/2. If site 2x is occupied and site 2x + 1 
empty, then the particle moves with probability 1 to site 2x + 1, i.e., 

T2x{t') = T 2x (t)T 2x+1 (t) 

(1) 

T 2x+ l{t') = r 2*(*) +T 2x +l{t) ~ T 2x {t)T 2x+1 {t) . 

These rules are applied in parallel to all pairs except the pair (L,l). In this pair representing 
the boundary (site L) and the origin (site 1) resp. particles are absorbed and injected 
according to the following stochastic rules. If site 1 was empty at time t then it remains 
so with probability 1 — a and becomes occupied with probability a at time t'. If site 1 
was occupied at time t then it remains occupied with probability 1. These two rules are 
independent of the occupation of site L. On the other hand, if site L was occupied at time 
t it remains so with probability 1 — (3 and becomes empty with probability f3. If site L 
was empty, it remains empty with probability 1. These two rules are independent of the 
occupation of site 1. This means that opposed to the models with sequential dynamics 
studied in refs. P,^,|5|,[ID|| simultaneous injection and absorption is allowed with probability 



a p. We have 

Ti(t') = 1 with probability Ti(t) + a(l — Ti(t)) 
Tl (f) = with probability (1 - a)(l - n(t)) 

(2) 

T L (t') = 1 with probability (1 - (3)r L (t) 

T L (t') = with probability 1 - (1 - P)r L (t) ■ 

In the second half step £ + 1/2— > t + 1 the pairing is shifted by one lattice unit such that 
the pairs are now (1,2), (3,4), ... [L — 1,L). Here rules (|I|) are applied in all these pairs, 
there is no injection and absorption in the second half time step.Q 



"Note that we reverse the order of the choice of pairs as compared to JTI. There the pairs were 
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In the mapping of ref. || this model is equivalent to a two-dimensional four- vertex model 
in thermal equilibrium with a defect line where other vertices, not belonging to the group 
defining the 6-vertex model or 8-vertex model, have non-vanishing Boltzmann weights. The 
two steps describing the motion of particles define the diagonal-to-diagonal transfer matrix 
T(a,fl) in the vertex model (see appendix). The pairing is chosen as in |5| but the hopping 
probabilities are different. 

For the Bethe ansatz solution of the p/p- model with a defect Q the conservation of 
the number of particles was crucial and we cannot repeat the calculation here, where the 
particle number is not conserved. However, since we are only interested in the steady state, 
we can construct the steady state explicitly for small lattices and then try to guess its general 
form for arbitrary length L. This method was succesfully applied in refs. [|J and || and 
led to exact expressions for the particle current and the density profile for arbitrary values 
of the injection and absorption rates. Only after guesswork produced the correct results, 
they were actually proven (see also |]IU||). It turns out that also here we can guess rules 
for the construction of the steady state. Instead of proving them we verified our conjecture 
for lattices of up to 14 sites. In the same way we guessed and verified expressions for the 
density profile (the one-point density correlator (r x ) in the steady state) eq. flTSP and the 
equal-time n-point density correlation function (|30D . The simple form of these correlation 
functions then allowed for a conjecture of the time-dependent two-point function ( T x T l T y ) 
(37|) - PH). (In this expression T* denotes the i-th power of the transfer matrix T(a,/3)). 
The mapping to the vertex model allows for an independent verification of this result. 

Among other things the time-dependent two-point function is of interest for the study 
of local dynamical scaling in the absence of translational invariance. Dynamical scaling in 
a 1+1 dimensional system with translational invariance both in space and time direction 
implies that the two-point function G(r,t) behaves under a global rescaling A of the space 



chosen as (1,2), (3,4), ... in the first half time step. 
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and time coordinates as 
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G(Xr, X z t) = \- 2x G{r,t) . (3) 

In this expression r denotes the distance in space direction, t is the distance in time direction, 
z is the dynamic critical exponent and x is the scaling dimension. From (|3[) follows that the 
correlation function has the form 

G(r,t) =r 2x/z §(^) (4) 

with the scaling function $ which is not determined by global dynamical scaling. By ex- 
tending the concept of global rescaling to local, space-time dependent rescaling, it has been 



shown that for the special case z = 2 the correlation function G(r, t) is of the form [0,0 



for 2 

G(r,t) = at~ x e'^t (5) 

with some constants a and b, i.e., $(r 2 /t) = aexp (— br 2 /t). The (connected) density cor- 
relation function in the probabilistic symmetric p/p-model without defect computed in |5j 
is indeed of this form with critical exponent x = 1/2. Since we study the steady state we 
have translational invariance in time direction, but due to the open boundary conditions 
translational invariance is broken in space direction. In sec. 5 we show that the form of 
( T x T^Ty ) for large L on the critical line a = (3 and in the scaling regime close to it resembles 
(|) with x = (sec. 5), i.e., one has z = 2, but there are additional pieces that arise from 
the breaking of translational invariance. 

The paper is organized as follows. In sec. 2 we present our conjectured rules for the 
construction of the steady state and exact expressions for the current and the density profile. 
In sec. 3 we study the limit of large L and derive the phase diagram. In sec. 4 we present 
expressions for the equal-time n-point density correlation function. They turn out to be 
reducible to a sum of one-point functions through associative fusion rules of the density 
operators. In particular we study the two-point function in the scaling regime. In sec. 5 
we compute the time-dependent two-point correlator. Again we put our emphasis on the 



vicinity to the phase transition line. In Sec. 6 we summarize our main results and discuss 
our results in the context of other exclusion models. In the appendix we discuss the mapping 
to a two-dimensional vertex model. 



II. CONSTRUCTION OF THE STEADY STATE 

Before we discuss the construction of the steady state we introduce some useful notations. 
In anticipation of the correspondence of the model to a vertex model discussed in refs. 
and in the appendix we denote a state of the system with N particles placed on sites 



x\, . . . , xn and holes everywhere else by \ x%, . . . ,Xn). The transfer matrix T(a, (3) (|A1|) 
acting on the space of states spanned by these vectors acts as time evolution operator and 
encodes the dynamics and the boundary conditions of the system as defined by eqs. (|]) 
and (0). The steady state is the (right) eigenvector with eigenvalue 1 of the transfer matrix 
T(a,j3) and we denote it by 

L 

1 1) = • • -, x n)\xi, ...,x n ). (6) 

N=0 {x} 

Here the iV-particle "wave function" $?n(xi, . . . , xn) is the unnormalized probability of 
finding the particular configuration \ xi, . . . ,Xn) of N particles in the steady state. We 
denote the state with no particles by | ) and the corresponding wave function by ^q. The 
summation runs over all states of N particles (0 < N < L) and all possible configurations 
{x} = {xi, . . . , xjv} and one has T(a, 1 } = | 1 ). The normalized probabilities are given 
by 

p N {xt, ...,x N ) = ^tv(xi, . . . , x N )/Z L (7) 

with 

Z L = • (8) 

^=0 {x} 

The transfer matrix T(a,(3) has a left eigenvector ( 1 1 with eigenvalue 1 given by 
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(i|= EEK-.^I . ( 9 ) 

N=0 {x} 

where ( x±, . . . , xjy | is the transposed vector to | xi, . . . , xn )• Defining a scalar product in 
the standard way (i.e., ( x 1: . . . ,xn \yi, . . . , x/m ) = 1 if the sets {x} and {y} are identical and 
else) one can write Zl as the scalar product Zl — ( 1 1 1 )■ 

Furthermore we denote the projection operator on particles on site x by t x : 



T x \xi,...,X N ) 



\xi,...,x N ) if x e {xi, . . . ,x N } 

(10) 

else 



The projector on holes is a x = 1 — r x . Expectation values ( t Xi . . . r Xk ) of the operators r x 
and their products in the steady state can conveniently written in the form 

(r Xl ...r Xk ) = (1\t Xi ... r Xk | 1 )/Z L . (11) 

Taking the scalar product with the left eigenvector ( 1 | and dividing by the normalization 
sum Zl is equivalent to a summation over all probabilities Pn(vi, • • • , Vn) with {x±, . . . , Xk} G 
{yi, . . . , i/n}- This is the definition of an expectation value in the steady state. 

The particle current j is a conserved quantity in the bulk since only the origin and the 
boundary act as a source or sink of particles. It is given by the correlator JT] (see ( |A6|) ) 

j = (T 2x a 2x+ i ) . (12) 

Now we discuss the construction of the steady state. In |l] we derived the important 
result that for the deterministic dynamics defined by (|l|) one has 

^-1^1 A) = (13) 

for 1 < x < L/2 and x < y < L/2 and any right eigenvector | A) of the transfer matrix. 
This simplifies the construction of the steady state considerably: If in a state \ xi, . . . ,xn) 
one of the X{ is odd, then it has a non- vanishing weight ^^(xi, . . . , x^) only if all even Xj 
with Xi < Xj < L are also contained in the set {xi, . . . ,xn}- 
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Using this it is easy to construct the steady state explicitly for small L. We discovered 
that the unnormalized probabilities (xi, . . . , x^r) in the chain with L + 2 sites can be 

constructed recursively out those of the chain with L sites according to the following rules: 



Rule 1: (0 < N < L, all {x}) 

¥^ +2 \ Xl + 2,x 2 + 2,...,x N + 2) = /3 2 (1 - a)*%\ Xl , x 2 , . . . , x N ) 

Rule 2: (0 < N < L, all {x}) 

¥^\2, xi, ... , xjv-i, L + 1, L + 2) = a\l - (3)^f{2, x u . . . , x N ^) 

Rule 3: (0 < N < L/2, {x u . . . , x L/2 } ^ {2, 4, 6, . . . , L - 2, L}) 

^Si 2) (2, Xl + 2, x 2 + 2, . . . , x N + 2) = a(3 2 ^ ] (x x , x 2 , . . . , x N ) ( 14 ) 

Rule 4: (L/2 < N < L, {x 1 , . . . , x L/2 } ^ {2, 4, 6, . . . , L - 2, L}) 

^N+i\ x i> x 2» ...,x n ,L + 2) = a 2 (3^\x 1} x 2 , ...,x N ) 

Rule 5: ({xi,...,x L/2 } = {2, 4, 6, . . . , L - 2, L}) 

*gg ) 1 (2, 4, . . . , L, L + 2) = a(3(a + /3)*S(2, 4, . . . , L) - (a/3) L / 2 + 3 



These rules together with (|T3| ) and the initial conditions 

v4 2) =/3 2 (l-a), ¥< 2) (2) = a/3, 2) = a\l - (3) (15) 

define recursively all quantities ^^ +2 \xi, . . . , x^) in the chain with L + 2 sites. 

Based on these rules we constructed the steady state up to L = 14 and verified that it has 



indeed eigenvalue 1 of T(a, f3). In a next step we computed the sum over all *&n(xi, . . . , xn) 
and concluded that the normalization Z L (R|) is given by 



L/2-1 



;i-/?k e 



fc=0 



a 



P)a 



L+l 



L/2-1 / \k 

(1 - a)p L+l 



a — P 

a L (l + L(l-a)) 



+ (af3) L / 2 

a^p 
a = P 



(16) 



This result was again checked explicitly up to L = 14. 

Going one step further we consider the average density (t x ) at site x defined by (PI) . 
We found the following exact expressions for the even and odd sublattices resp.: 

2x-l / n\k L 



Zl{t 2x ) 



E (rV + ^-^'E f-V + M) L+1 



k=0 



a 



k=2x 



a 



a L+1 (1 + L{1 - a)) + 2xa L {l - af 



(17) 



Z l (t 2x -i) 



2x-S 

E 

fc=0 



i - pf a L E ( - ] + (i - '*K :;Z 



L+2-2x n2x-2 



a 



L+l, 



1 - a) + (2x- l)a L (l - a) 



On the r.h.s. of (|17]) the upper expressions are valid for a ^ P, while the lower expressions 
are valid for a = p. Performing the summations we arrive at the main result of this section 

1 _ (H)2x 



>2x) 



a + (1 — a) 



a + (1 — a) 



i _ l = a_(P\L+ 

2x 



1 + L(1 - a) 



- a^P 

a = P 



(t 2x -i ) 



a 



i _ l-a { /3\2x-l 

1 _ l=a(§_\L+l 
1 1-/3 I a ) 



(18) 



2x 



+ 



ail — a] 



a^P 
a = P 



l + L(l-a) l + L(l-a) 
The anisotropy between the even and odd sublattices is a consequence of the parallel 

updating mechanism JTJ and related to the net particle current flT2| ) for which we found 

(again, by direct evaluation on small lattices and guessing the general form for arbitrary L) 
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J = s 



a 



1 _ i-a(£\L+i a "r P 



a(l — a) 



(19) 



a — 



a = (3 



1 + L(1 - a) 

The fact that this quantitity is independent on x as it should be is an additional non-trivial 
check for our conjectures. 

Finally note that by the definition of the model there is the particle-hole symmetry ( |A5[ ): 
Changing particles into holes, reflecting site x into site L + 1 — x and exchanging a and j3 
leaves the system invariant. All our results are indeed invariant under this operation. 



III. THE PHASE DIAGRAM 



From eqs. (18) one realizes that the system changes its behaviour if a = /3. The average 
densities on the even and odd sublattices 



L/2 



L/2 



p m = 7 E(^). p (odd) = 7 E(^-i) 

L x=l L x=l 



(20) 



have a discontinuity in the thermodynamic limit L — ► oo at a = (3 ^ 0, 1.Q One finds from 
(0) 



P 



(even) 



(odd) 



a a < j3 

1 a > (3 

a < (3 

I- (3 a> (3 



(21) 



tlf a or (3 is either or 1 the system is trivial in the sense that the density profile is exactly constant 
(even in finite systems) and all correlation functions can be obtained without any calculation (see 
appendix). Therefore we exclude these cases from our discussion. 
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On the phase transition line a — (3 one obtains p( even ) = (14- a)/2 and p( odd ) = (1 — a)/2 
resp. For a < (3 (more particles are absorped than injected) the system is in a low density 
phase with total average density p—\-, (p( even *> + p( oa!a! )) = a/2 < |, while for a > (3 it is in 
a high density phase with p = 1 — (3 /2> \ (Fig. 1). 




a 



FIG. 1. Phase diagram of the model in the a — (3 plane. Region A is the low density phase and 
region B the high density phase. The phases are separated by the curve a = (3. 

In the thermodynamic limit L —>■ 00 the current j ([19|) is given by 

j = min (a, (3) . (22) 

There is no discontinuity at a = (3 in the current, but its first derivatives w.r.t. a and (3 are 
discontinous. The discontinuity of p at the phase transition line a = (3 and eq. (0) remind 
us of the s/o- model 0,^,0]. In this model in the region a, (3 < 1/2 the phase diagram shows 
a low density phase Ai and a high density phase Bi separated by a phase transition line at 
a = (3 [[J. Also in this model the density and the first derivatives of the current w.r.t. the 
injection and absorption rates a and (3 have a discontinuity at the phase transition line. 

In terms of the sublattice densities the current j is given by j = p^ even ^ —p(° dd ) for all a, (3. 
In terms of the total average density p the current satisfies j = 2p if p < 1/2 (low density 
phase) and j = 2(1 — p) ifp> 1/2 (high density phase). These are the same relations as in 
the p/p-model in the respective phases [l|]. 

10 



Now we turn to a discussion of the density profile. We first study the case a < (3 and 
L — > oo. Defining the decay length £ by 



T^ln^ (23) 

one obtains from ( |T8| ) the density profile up to corrections of order exp (— L/£) 

(r 2x ) = a + (l- / 9) e -( i + 1 - 2 ^ 

(24) 

= (1 - /5) e -(^ 2 - 2 ^ . 

The profile decays exponentially with increasing distance from the boundary to its respective 
bulk values p^uf 1 = a an d Pbuik = 0- This the low density phase of the system. 

In the high density phase a > (3 which is related to the low density phase by the particle- 
hole symmetry the profile is given by 

(t 2x ) = l-(l-a)e- 2 ^ 

(25) 

The bulk densities are p[^ n) = 1 and p[°J° = 1-/3. 

On approaching the phase transition line a = f3 the decay length £ diverges. On the line 
the profile is linear and up to corrections of order L _1 given by 

2x 

(26) 



(r 2x ) = a + (1 - a) : L 



2x — 1 

(T2x-i) = (1-a) L ■ 

An explanation for the shape of the profile in the two phases and on the phase transition 
line will be given in the next section. 



IV. EQUAL-TIME CORRELATION FUNCTIONS 



Having found exact expressions for the current and the density profile we proceed calcu- 
lating the n-point equal-time density correlation function (t Xi . . . r Xn ) in the steady state. 
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Examining the two-point function for small L we found the following exact relations for 

1 < x < y < L/2 

(r 2x r 2y ) = (r 2x ) - a + a(r 2y ) 
(r 2x r 2 y-i) = (1 - (3){{t 2x ) -a) + a{r 2y -i) 

(27) 

(r 2x+ iT 2y ) = (r 2x+ i) 

(T 2x -lT 2 y-l) = (1-P){t2x-i) ■ 

The third of these equations is a simple consequence of (|13"D which says that whenever there 
is a particle on an odd lattice site then all even lattice sites to its right must be occupied 
as well. The important result is that the two-point is completely determined by the one- 
point function and some constants! Going further we made the surprising observation that 
the n-point function can also be expressed in terms of one-point functions by repeatedly 
fusing products of operators r x r y according to the fusing rules that are defined by (|27|) by 
omitting the averaging. This fusion can be performed in arbitrary order until one reaches 
the one-point level. 

The fusion rules implied by eqs. fl2"7p can be simplified by using operators rj x defined by 

T 2 X — Ct T 2x _i 
V2 X = -j , = i a ( 28 ) 

1 — a 1 — p 

instead of using the density operators r x . In the bulk of the high density region both ( r] 2x ) 
and ( 7722-1 ) take the value 1, while in the bulk of the low density region both average values 
are 0. Expressing all r x in terms of the r] x the correlation functions ( p7[) become 

(*Mea) = (Vxi) (X 2 >X 1 ). (29) 

Fusion of n operators rj Xl . . . T] Xn gives i] x . with X{ = min {x±, . . . , x n }. So the n-point corre- 
lation function is 

(Vxi ■■■Vx n ) = (Vxi) (Xi = min{xi, . . . ,x n }) . (30) 

This is the main result of this section. 
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The form of the two-point function fl29|) can be understood by considering the steady 
state as composed of "constituent profiles" with a region of constant low density up to some 
point Xo in the chain followed by a high density region beyond this "domain wall". Such 
an assumption explains why the correlator does not depend on X2- In the low density 
region of density a on the even sublattice and on the odd sublattice the operator rj xi has 
vanishing expectation value and therefore the whole expression (t] Xi t] X2 ) is zero if x\ is in 
this region, independent of r] X2 . If, however, x\ is in a region of high density, then, according 
to our assumption, also x<i > x\ must be in region of high density. Thus, r\ Xl r\ X2 again does 
not depend on xi and takes the value 1. We conclude that the product f] Xl rj X2 is either or 
1, depending on whether x\ is in a region of low or high density. This leads to the expression 
(f29p for the expectation value of this product. 

The average value ( r\ x ) itself contains the information about the position x$ of the 
domain wall. In the low density phase the density profile decays exponentially from above 
to its bulk value with increasing distance from the boundary. This means that the probability 
of finding the domain wall also decreases exponentially with the same decay length £ with the 
distance from the boundary. The domain wall is caused by particles hitting the boundary 
where they get stuck with probability 1 — (3 and then cause other incoming particles to 
pile up and create a region of high density (fig. 4 in the appendix). On the other hand, in 
the high density phase where the rate of injection is higher than the absorption rate the 
situation is reversed. Here the probability of finding the domain wall decreases exponentially 
with the distance from the origin. This can alternatively be explained either in terms of 
holes through the particle-hole symmetry or in terms of particles being piled up from the 
boundary over the whole system up to a point close to the origin. On the phase transition 
line injection and absorption are in balance and the probability of finding the domain wall 
is space-independent and the density profile is a linear superposition of the assumed step 
function constituent profiles. This leads to the observed linearly increasing average density 

(0)- 

We conclude this section by studying the two-point function in more detail in the limit 
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L — > oo. We define the equal-time connected two-point function by 

G c ( Xl ,x 2 ;t = 0) = ( VxiVx2 )-(v* l )(v* 2 )= { Tx ^ X 2 ^~ ^ } ■ ( 31 ) 

where m = 2, n = if both x± and x 2 are even, n = 2, m = if both x± and x 2 are odd and 
m = n = 1 else. 

In what follows we restrict ourselves to the case where X\ and x 2 are both odd, the mixed 
correlators can be computed analogously. From ( [29D one obtains G{x\,x 2 \^S) = (r) xi )(1 — 
(Vx 2 )) an d inserting the expressions for the density profile (p~3|) - (|T5| ) one obtains with 
x 2 = Xi + 2r (r > 0) 

' A{x 2 )e- 2r/ t a<(3 
G c (x u x 2 ;0) = I A( Xl )e~ 2r ^ a>(3 • (32) 

^(l-i)-fi a = /9 
The amplitudes of the exponential decay are given by 



A(x) = e~ R /t (l - e- R /« 



(33) 



where i? = L + 1 — x measures the distance of site a; from the boundary. 

The decay length £ is identical with correlation length of the connected two-point func- 
tion. On the phase transition line the correlation function is constant for relative distances 
2r <C L. Its amplitude depends on the position x\ in the bulk. A similar form of the 
connected equal-time correlator was found in the p/p- model with a defect 

V. TIME-DEPENDENT CORRELATION FUNCTIONS 

In this section we study the time-dependent two-point correlation function in the steady 
state 

G(xi, x 2 ; t) = ( rj Xl T l rj X2 ) (34) 
14 



where T* denotes the t-th power of the transfer matrix T.f] We define the direction of the 
time evolution formally by T~ t T x (to)T t = r x (t + t), thus G(x 1 , x 2 ; t) = ( r] Xl (t + t)rj X2 (t ) ). 
The connected two-point function is defined by G c (xi,X2] t) = G(xi, x<i\ t) — ( r) Xl ) ( r\ X2 ). 

The standard way of computing the correlation function fl3"lD would be the insertion of 
a complete set of eigenstates of T, evaluating the matrix elements a^xi) = ( 1 \f] Xl \ A*. ) and 
SfcO^) = ( Afc \r] X2 \ 1 ) and summing over afca&A^.. Since we do not know the eigenstates and 
eigenvalues we take the alternative route using the commutator of [i] x , T*]. Since ( 1 | is a 
left eigenvector of T with eigenvalue 1 one has ( [t) x , T l }r] X2 ) = ( 1 K^T* — T t r] xl )r] X2 \ 1 ) = 
G c (xx, xi\ t). From this one obtains G{x\, x%; t). 

First we note that from the commutation relations ( |^7j ) of r x with T one obtains 

T2x-1 T = T (1 — <J2x-2<3'2x-l) T 2x T 2x+l 

(35) 

T 2xT = T (1 — 0"2x-20"2a;-l(l — T 2x T 2x+l)) ■ 

It is obvious that evaluating t x T* is not an easy task. By iterating relations (|35|) t times 
not only the number of terms in the products on the r.h.s. but also the total number of 
such multi-point correlators increases extremely fast with t. It is only the simplicity of the 
multi-point correlators (see eqs. (p9|) and Q30D) that makes this approach promising. We 
restrict our discussion again to both x\ = 2yi — 1 and X2 = 2y2 — 1 odd. 
By iterating ( p5[ ) t times one finds that T2 yi -\T t T2y 2 -i is of the form 

T~2yi-lT t T2y 2 -l = T* { 1 — (o"2yi-2t cr 2yi-2t+l • • • ) — ( cr 2y 1 -2t+2°~2yi -2t+3 ■ ■ ■ )~ 

( • • • ) — {o'2y 1 ~2t+2kO'2y 1 -2t+2k+l ■ ■ ■ ) ~ ( • • • )~ (^) 
0" 2y 1 -4+2tO~2yi -3+2t } T~2 yi -2+2tT~2 yi - l+2tT~2y 2 - 1 

where the dots denote some complicated sums of products of operators r y acting on sites y 
between 2yi — 2t and 2yi — 1 + 2t. a y = 1 — r y denotes the projector on holes and in order to 



^Note that t = 1 corresponds to a distance of two lattice units in time direction in the underlying 
vertex model (see appendix). 
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avoid boundary effects one has to choose t < y x — 1. We first discuss the correlation function 
outside the light cone, then on the edges of the light cone and finally in its interior. 



A. Correlation function outside the light cone 

We want to evaluate (T2 yi - X T t T2y 2 -i ) with 2y x — 1 > 2y 2 — 1 + 2t. Recalling the fusion 
rule ( T2x+iT~2y ) = ( t 2x + x ) for y > x (|27|) one obtains ( T 2x +iO'2y ) = for y > x. Since 



the fusion procedure is associative all terms on the r.h.s. of (p6|) vanish when contracted 
with T 2y2 -i except ( T 2y2 _iT 2yi _ 2+2t T 22/1 _ 1+2 i ) = ( T2y 1 -i-2tT2 yi -2+2tT2y 1 -i+2tr angle . Using 
also ( T2 X -iT 2 y-i ) = (1 - P){t 2x -i) for y > x one gets ( T 2m _i_ 2 tT 2yi _ 2+2t T 2ai _ 1+2 i ) = (1 - 
P){ T2yi-i-2t )• Therefore we obtain 

G(xi,x 2 ;t) = (r] X2 ) (x x , x 2 odd, x x > x 2 + 2t) . (37) 

Now we study the correlator ( T 2 j, 1 _iT t T2j /2 „i ) with 2y x — 1 < 2y 2 — 3 — 2t. Here the 
fusion of r 2yi -i + 2t in the r.h.s. of (|36|) with r 2j/2 _i yields (1 — /5)r 2j , 1 _ 1+2t and by taking the 
average value one obtains ( T2j /1 -iT*t 2 j /2 _i ) = (1 — 0){ T2 Vl -\T t ) = (1 — /3)(r 2yi _i ). We find 

G(xi,x 2 ;t) = (rj Xl ) (x x , x 2 odd, x x < x 2 - 2 - 2t) . (38) 

Eqs. (^) and (|38| ) are no surprise. The area defined by (137|) and (|38D is the exterior of 
the forward light cone of the particle at site x 2 . If x x and x 2 are chosen in this way and 
both are in a region of uniform density (either in the bulk of the high density phase or in 
the bulk of the low density phase) one has 7] Xl = r] X2 = 1 or and therefore the connected 
correlation function G c (x x , x 2 ; t) is time-independent and as one would expect. In the 
boundary region of the low density phase where particles pile up and lead to a non-uniform 
density profile (or near the origin in the high density phase) it is still time-independent as 
it must be outside the light cone, but non-zero (see ([32])). This is due to the hard-core 
repulsion of the particles which behave as an incompressible liquid. 
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B. Correlation function on the edges of the light cone 

On the right edge of the light cone of the particle at site X2 defined by x\ = X2 — 
2 + 2t we can repeat the considerations that led to (|37j): All the pieces on the r.h.s. of 
(H) containing T 2y2 -i<72 yi -2t+2k = T 2yi+ i-2t^2 yi -2t+2k with k > 1 vanish as a result of the 
fusion rules and only the first two pieces in the sum remain. Although the term containing 
°"2j, 1 -2<o"2j/ 1 +i-2t • • • T2yi+i-2t does not vanish due to fusion with r 2yi+ i^2t it is nevertheless 
since by definition o 2yx +\-2tT2 y ^+\-2t = 0. Therefore 

G(x 2 -2 + 2t,x 2 ;t) = (r] X2 ) (x 2 odd) . (39) 

Consequently the connected correlation function on the odd sublattice vanishes also on the 
forward edge of the light cone if the two points are in a region of uniform density. This is 
a result of the asymmetry of the model: if the system is in a region of uniform low density 
p < 1/2 the odd sublattice is empty and the vanishing of the correlator is trivial. In a 
region of uniform high density the even sublattice is completely occupied and particles on 
the odd sublattice effectively move only to the left (fig. 4 in the appendix) and are therefore 
uncorrelated to particles on the right edge of their (forward) light cone. 

Due to the deterministic nature of the dynamics the particles on the odd sublattice 
move with the velocity of light, i.e. two lattice units per full time step as long as they 
are in region of uniform high density. Thus we expect a singularity of the correlation 
function on the left edge of the light cone defined by x\ = x 2 — 2t: Indeed, choosing 
2?/2 — 1 = 2yi — 1 + 2t does not change the r.h.s. of (|36|) since Tf yi _ 1+2t = r 2yi ^i +2t and 
therefore (T2 yi -iT t T 2 y 1 -i+2t ) — ( r 2?/i-i^') = ( r 2 yi -i)- For the correlator (j34|) we obtain 

G(x 2 -2t,x 2 ;t) = (l-P)- 1 (r ]x2 _ 2t ) (x 2 odd) . (40) 

Here the connected correlation function in a region of uniform high density does not vanish. 
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C. Correlation function inside the light cone 



First we note that in a region of uniform high density the result is again trivial. In such 
a region particles on the odd sublattice are found everywhere with equal probability (the 
equal-time connected two-point function is 0) and since they move with the velocity of light 
the time-dependent connected two-point function does also vanish. 

If the profile is not uniform the calculation inside the light cone is non-trivial. With 
x\ = 2yi — 1 as above and x 2 increasing beyond 2yi + 2 — 2t more and more contributions 
from the r.h.s. of fl3"E| ) are non-zero. We evaluated G(xi, x 2 ; t) for t = 1, 2, 3 inside the light 
cone on the computer (using the software system Mathematica IfLljl ) by calculating the exact 
form of (|36|) and then implementing the fusion rules (27) on the multi-point correlators 



on the r.h.s. of (|36|) . First we noticed that in a region of uniform high (low) density (all 
(rj x )=l(0)) one obtains G(xx,x 2 ;t) = 1(0) and therefore G c (xi, x 2 ; t) = as it should be. 
This observation is indeed a highly non-trivial test of the conjectured fusion rules (|27|) on 



which our calculation is based: Since eqs. ( p?[) are supposed to be exact, the result of the 



calculation of the time-dependent correlator must also be exactly 1 (0) if all rj x involved are 
set to 1 (0). Any other result would have shown that the fusion rules do not hold. Secondly 
we observed that by taking a = 1 — (3 the exact general form of the correlator becomes fairly 
obvious for arbitrary values of x%, x 2 and t with x\ inside the light cone of x 2 . We found by 
generalizing our result from t = 1, 2, 3 to arbitrary t 



G(x,x + 2y;t) = ^ j ^/3 2 ^(1 - ft 



))}+ (41) 

This the main result of this section, valid for t < (x — l)/2 and —t + 2<y<t—l. The 
first restriction is due to boundary effects, the second defines the interior of the light cone. 
In order to check this result we explicitly calculated G(x, x + 2y; t = 4) on the computer 



for arbitrary a and ft using ([36]) and the fusion rules and found it in exact agreement with 
our conjecture when setting a — 1 — ft. As a second, independent test we set r) x = 1(0) 
corresponding to the bulk value in the high density region (low density region) and indeed 
obtained G c (xi, x 2 ] 4) = for arbitrary a and ft. 

The choice a = 1 — ft is not too restrictive as far as the physics is concerned: since 
this curve runs across the phase diagram it covers both the high density phase and the low 
density phase and crosses the phase transition line at a = ft = 1/2. In what follows we 
study G(xx,X2',t) in the low density phase along the curve (3 = 1 — a > 1/2 and on the 
phase transition line at ft — 1/2. 

In the low density phase we focus on the boundary region with a non-uniform density 
profile. For ft > 1/2, a = 1 — ft the expression fl23|) for the density profile for large L gives 
( tfcs-i ) = (V2 X ) = ((l- ft)/ft) L+2 - 2x and therefore 

/, a \ L+l-x+2t-2k 

ft( Vx+3 . 2t+2k ) + (1 - ft)(Vx+4-2t+2k) = \^f) • (42) 



Inserting this into (|4l| ) and introducing the incomplete /3-function 

t+y-2 / \ 

E r: 2 )p 2t - 2 - k (l-P) k = I fi (t-y,t + y-l) 



k=0 \ k 



(43) 



1 - h.pit + y,t-y)- r^jft^il - ft) t+y - 1 



the correlation function ([41~D can conveniently be rewritten 



G(x,x + 2y,t) = (±fi) [Ii-f,(t-y,t + y) + ( 1 ^) h^(t + y,t-y)J . (44) 

For large times t (such that \y\/t <C 1) the incomplete /5-function has the asymptotic 
form 



h-p{t + y,t-y) 



(45) 



with 
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C 1 = - In (4/5(1 - /3)) f C 1 = - In ^ - (46) 

and the probability integral P(u) = 1/V2n /"^ exp (—t 2 /2)dt. In terms of £ t the inequality 
/3 > 1/2 in the upper expression of the r.h.s. of ([45]) has to be understood as 1 £t^$t. In 
the lower expression we assume 1 <C t^$£f. Note that the two length scales £t and £ r are not 
independent quantities but related through £ f -1 = lncosh 2 (^~ 1 /2). As (3 approaches 1/2, £ f 
and £ r diverge and are asymptotically related through £ f 4£ 2 . 

We define £ = £ r and r = |y| = l/2|x2 — x\\ and insert (^5j) into (fH|). This gives the 
scaling form of the time-dependent correlation function in the scaling region of large 2£i$i 1 / 2 



G(xi,x a ;t) = e-Me-'Kj^e-VWe-'*'* (47) 

V 7rt 



where 

L + 1 — X2 if 2y = 22 — a?i > 
i2 = < (48) 

L + 1 — 2i if 2y = x 2 — x\ < 

measures the distances of x 2 or x\ from the boundary, depending on the sign of x 2 — X\. 
G(xi,X2',t) is invariant under the scaling transformation 

R — > Ai?, r -> Ar, ^ -> A£ 

(49) 

t -> A 2 t 

This is of the form corresponding to dynamical scaling with a dynamic critical exponent 
2 = 2 and critical exponent x = 0. 

If £ increases beyond the crossover length scale t 1//2 the correlation function changes its 
form. At the critical point f3 = 1/2 we have (rj x ) = x/L and up to corrections of order 1/L 
the exact expression ( 4~I|) for the correlation function gives 



x + 2y 2t + 2y T . 
G(x,x + 2y;t) = — ^ (t - y,i + y - 1) + 
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Using 

t+y-2 



2 £ tl 



(3 2t-2- k{l _ (3)k = 2(i-p)(2t-2)If ) (t-y,t + y-l)- 

(51) 



2(t + y-l)!^-_ 2 i j^(l-/3)^- 1 



and the expansion ( f45| ) of Ii(t + y,t — y) we obtain with 2y = x 2 — X\ and u = y/y/i 



2 



2 



(52) 

Xl+X2 ^ fe-" 2 + v ^n(l-2P(n)r 



2L v^FL 

For finite distances y, t one has (up to corrections of order L~ l which we neglect) x\ = x 2 = x 
and the correlation function at the critical point is space and time independent with an 
amplitude x/L depending on the relative position of x in the bulk. For large times, t oc L 
and y/t <C 1, the correlation function gets an contribution order L~ x l 2 . G{x,\^x 2 \t) is 



invariant under the scale transformations P9| ) with the length scale £ replaced by the size 
of the system L. 

This result has a simple interpretation in terms of the constituent profiles discussed in 
the preceding section. For simplicity we consider y = 0. The operator ^T*^ gives 1 if 
x is in a region of high density both at times to and to + t. The probability that x is in 
the high density region at time to is x/L. This accounts for the constant x/L in ( |52"|) If we 
assume that the domain wall separating the low density region from the high density region 
performs a random walk of two lattice units per time step starting from its position xq at 



time to then the resulting expectation value {r/ x T rj x ) will indeed be of the form (52). 



VI. SUMMARY AND COMPARISON WITH OTHER EXCLUSION MODELS 

Let us summarize our main results. We obtained recursion rules (|nj) which allow the 
construction of the steady state of the model defined by ([]]) and (0). Moreover we found for 
arbitrary values a and j3 of the injection and absorption rates exact steady state expressions 
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for the density-profile (|T8|), the current (|19D, the equal-time n-point density correlation 
function ( |30| ) and the time dependent two-point function eqs. ( |37| ) - ([4TD f\ 
We made the following observations: 

(a) The phase diagram (fig. 1) shows two phases. In the low density phase (a < (3) the 
average density is p = a/2 (in the limit L — > oo) and in the high density phase one has 
p = 1 — (3/2. On the phase transition line a = (3 the average density is p = 1/2. 

(b) There is a shift in the average densities between the even sublattice and the odd sublat- 
tice (see eqs. fl24|) - (|26)). This shift is the current j. In terms of the total average density 



p one finds j = 2p in the low density phase and j = 2(1 — p) in the high density phase. On 
the phase transition line (p = 1/2) the current is j = a (i.e. j ' ^ 1 = 2p = 2(1 — p)). 

(c) In the low density phase the density profile decays on both sublattices to its respective 
bulk value exponentially with increasing distance from the boundary while in the high den- 
sity phase it increases exponentially to its respective bulk value with increasing distance 
from the origin. At the phase transition line a = (3 the length scale £ ( |23| ) associated with 
the exponential shape of the profile diverges and the profile increases linearly on both sub- 
lattices. The average density and the first derivatives of the current w.r.t. a and (3 have a 
discontinuity at the phase transition line (in the thermodynamic limit L — > oo). 

(d) The decay length £ is identical with the correlation length of the connected two-point 
density correlation function. Outside the light cone this correlator is of the scaling form 
G(xi,X2',t = 0) = Ar K exp (— r/£) with exponent k = 0. The amplitude A is space- 
dependent as a result of breaking of translational invariance fl32|) and (0). (e) The time- 
dependent two-point correlation function G(xi,X2',t) inside the light cone near the critical 
line ( f47D has a form compatible with dynamical scaling with dynamic critical exponent z = 2 
and scaling dimension x = (see the transformations (f49D). It contains the exponential 



§The expression ( f4l| ) for the time-dependent two-point function inside the light cone was found 
only for a = 1 — (3. 
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exp (— r 2 /t) characteristic for local dynamical scaling, but the amplitude is again space- 
dependent due to breaking of translational invariance. The correlation function changes its 
form when the correlation length increases beyond a crossover length scale of order t 1 / 2 . On 
the critical line it becomes a space-dependent constant up to corrections of order L _1 if t is 
large but finite, and up to corrections of order L -1 / 2 for times of order L. 
(f) From an analysis of the two-point correlation function we found that the density pro- 
file can be considered as a superposition of step-function type profiles with average den- 
sity p( et,en ) — a anc i p( odd ) — o up to some point xq and average density = 1 and 
P2 = 1 — (3 beyond this point up to the boundary. These densities are the sublattice 
densities in the low density phase and high density phase respectively. The probability of 
finding the "domain wall" at site x separating the two regions of high density and low 
density decreases exponentially with increasing distance from the boundary (origin) in the 
low density (high density) phase. On the phase transition line this probability is space inde- 
pendent and the domain wall can be found everywhere with the same probability. Studying 
the time-dependent correlation function suggests that it performs a random walk around its 
position x at time t . 

We conclude our discussion of the deterministic p/o- model with a brief comparison with 
other exclusion models and some conjectures for probabilistic exclusion models with parallel 
dynamics. 

In the deterministic p/p-model with a defect we found a phase diagram showing a low 
density phase, a coexistence phase where a low density density region coexists with a high 
density region and a high density phase. The relation between the current and the sublattice 
densities and the total density discussed in (b) for the low and high density phases is identical 
with that in the p/p-model in the respective phases. This is a consequence of the bulk 
dynamics which are identical in both cases and could have been guessed. 

It is interesting to observe that there is also a correspondence between the coexistence 
phase of the p/p-model and the phase transition line here: In the p/p-model with defect 
strength 1— q the quantity q is the hopping probability at a single link of the ring, say between 
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sites L and 1 and therefore corresponds to an absorption of particles at site L and injection 
of particles at site 1 with rate q. This injection and absorption is correlated because of the 
particle number conservation. In the coexistence phase the current is density-independent 
and one has j = q. In the p/o-model on the phase transition line discussed here particles are 
also injected and absorbed with the same rate a, but uncorrected. The current is density- 
independent, j = a as in the p/p- model. The profile is build by constituent profiles with a 
region of low density up to some xq and a high density region beyond that point, xq can be 
anywhere with same probability. This suggests that also the profile in the coexistence phase 
of the p/p- model with defect is build by such constituent profiles, with the distinction that 
there Xq cannot be anywhere (because of particle number conservation) but the probability 
p(xq) of finding the domain wall is centered around some point Rq. We believe that similar 
phases occur also in probabilistic p/p- models with a defect and in probabilistic p/o-models. 
In such models particles do not always move if the neighbouring site was empty as here but 
can stay with some non-zero probability even in the bulk. For mixed models (defect and 
some uncorrelated injection and absorption) we expect correspondingly a softening of the 
distribution p(x ). 

All the features summarized under (c) (except the implied anisotropy between the even 
and odd sublattices) are in common with the phase transition from the low density phase 
Aj to the high density phase Bi in the s/o-model |J. This suggests that also the correlation 
functions are qualitatively the same (in the thermodynamic limit L — > oo) near the phase 
transition line. In the s/o-model however, the transition line a = (3 extends only up to 
a = (3 = 1/2 as opposed to a = (5 = 1 here. Correspondingly, in our model there is 
no maximal current phase with a power behaviour of the density profile and no phases 
corresponding to the phases A n or B n (for a > 1/2 or (3 > 1/2) of where the shape of 
the density profile is determined by a product of a power law behaviour with an exponential 
decay. 

We believe that phase transitions to such phases cannot occur in our model because of 
the deterministic nature of the dynamics. In || we argued that these phases result from 
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an "overfeeding" of the system with particles: The system reaches its maximal transport 
capacity at density 1/2. In order to obtain average density 1/2 at the origin particles have 
to be injected with rate 1/2. If particles are injected at a higher rate they block each other 
rather than moving into the bulk and cause a phase transition. Here this cannot happen. 
The system reaches its maximal transport capacity also at average density 1/2 but here this 
corresponds to a completely filled even sublattice and an empty odd sublattice. An average 
density of 1 on the even sublattice at the origin can only occur if particles are injected at 
rate 1. Thus there can be no "overfeeding". The deterministic hopping rules imply that 
particles injected at the origin move away with velocity of light, therefore there can be no 
mutual blockage near the origin. (Similar arguments can be used for a discussion of the 
dependence of the phase transitions on the absorption rate (3 by exchanging particles with 
holes and studying the injection of holes at the boundary). 

This discussion naturally leads to the conjecture that the probabilistic p/o-model (which 
has not yet been studied) will have a phase diagram similar to that of the s/o- model with 
a phase transition line at a = (3 up to some value «o where the profile is constant on both 
sublattices and additional phase transition lines at the values of a and (3 corresponding to 
an overfeeding of particles at the origin and holes at the boundary respectively. 
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APPENDIX A: MAPPING TO A TWO-DIMENSIONAL VERTEX MODEL 

Following the idea of ref. || we want to show how the exclusion process defined by 
relations ([I]) and (0) is related to a two-dimensional vertex model. The discussion is partly 
similarly to that in [|1[] which we repeat for the convenience of the reader unfamiliar with 
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this mapping. The mapping of the boundary conditions to vertices in the vertex model is 
different from 

Consider a 4-vertex model on a diagonal square lattice defined as follows: Place an up- 
or down-pointing arrow on each link of the lattice and assign a non-zero Boltzmann weight 
to each of the vertices shown in figure 1. (All other configurations of arrows around an 
intersection of two lines, i.e., all other vertices, are forbidden in the bulk.) The partition 
function is the sum of the products of Boltzmann weights of a lattice configuration taken 
over all allowed configurations. 




"1 «2 <J2 ^2. 

FIG. 2. Allowed bulk vertex configurations in the four-vertex model. Tip-pointing arrows corre- 
spond to particles, down-pointing arrows represent vacant sites. In the dynamical interpretation 
of the model the Boltzmann weights give the transition probability of the state represented by the 
pair of arrows below the vertex to that above the vertex. 

In the transfer matrix formalism up- and down-pointing arrows in each row of a diagonal 
square lattice built by M of these vertices represent the state of the system at some given 
time t. Corresponding to the M vertices there are L = 2M sites in each row represented 
by the links of the diagonal lattice. The configuration of arrows in the next row above 
(represented by the upper arrows of the same vertices) then corresponds to the state of 
the system at an intermediate time t' = t + 1/2, and the configuration after a full time 
step t" — t + 1 corresponds to the arrangement of arrows two rows above. Therefore each 
vertex represents a local transition from the state given by the lower two arrows of a vertex 
representing the configuration on sites j and j + 1 at time t to the state defined by the 
upper two arrows representing the configuration at sites j and j + 1 at time t + 1/2. The 
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correspondence of the vertex language to the particle picture used in the introduction can 
be understood by considering up-pointing arrows as particles occupying the respective sites 
of the chain while down-pointing arrows represent vacant sites, i.e., holes. 

The diagonal-to-diagonal transfer matrix T acting on a chain of L sites (L even) of the 
vertex model with vertex weights cl\ , . . . , C\ as shown in fig. 2 is then defined by P|,pT5[ 



L/2 L/2 

T =n T 2 J -i-n T 2i= T ° ddTe 

i=l 3=1 



(Al) 



The matrices Tj act nontrivially on sites j and j + 1 in the chain, on all other sites they 
act as unit operator. All matrices Tj and Ty with \j — j'\ ^ 1 commute. (The difference 
j — j' is understood to be mod L). For an explicit representation of the transfer matrix 
we choose a spin- 1/2 tensor basis where the Pauli- matrix crj acting on site j of the chain is 
diagonal and spin down at site j represents a particle (up-pointing arrow) and spin up a hole 



(down-pointing arrow). In this basis Tj 



of) is the projection operator on particles 



on site j, Oj = |(1 + crj) is the projector on holes and sf = \{p X j ± ^0"J) (a x ' y ' z being the 
Pauli matrices) create (sj) and annihilate (s+) particles respectively. 

The bulk dynamics of our model is encoded in the transfer matrix by choosing the vertex 
weights as 



at = a.2 



In the bulk this leads to 



1 + Si S j+1 - T j( J j+1 



a = i 



1 \ o o 



110 

oooo 

1 



(A2) 



(A3) 



' 3,3- 



In the particle language the matrices Tj describe the local transition probabilities of particles 
moving from site j to site j + 1 represented by the corresponding vertices. If sites j and 
j + 1 are both empty or occupied, they remain as they are under the action of Tj. The same 
holds for a hole on site j and a particle on site j + 1, corresponding to the diagonal elements 
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of Tj, representing vertices ai, a 2 and c\. If there is a particle on site j and a hole on site 
j + 1, the particle will move with probability one to site j + 1. This accounts for vertex b 2 . 

As discussed in the introduction we assume open boundary conditions with injection of 
particles on site 1 and absorption of particles on site L. This allows for the additional vertices 
shown in fig. 3 together with vertex weights corresponding to the respective probabilities of 
creating and annihilating particles. 



FIG. 3. Additional vertex configurations allowed at the boundary and their Boltzmann weights. 
The left arrows of these vertices describe the particle configuration at the boundary site L of the 
system while the right arrows define the particle configurations at the origin (site 1). 

In a two-dimensional lattice (fig. 4) we consider the half-vertices at the left boundary as 
the right arms of the vertices shown in fig. 2 and fig. 3 and the half-vertices at the right 
boundary as their left arms. Thus the left arrows define the particle configuration on site L 
and the right arrows are considered as site 1. Vertices a±, a 2 and b 2 have a different weight 
at the boundary: a[ — 1 — j3, a 2 — 1 — a, b 2 = af3. Note that vertex b 2 at the boundary 
describes simultaneous absorption of a particle at site L and creation of a particle at site 1. 




a 
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1 2 3 4 5 6 7 8 9 10 11 12 
FIG. 4. Configuration of particles (up-pointing arrows) on a lattice of length L = 12 in space 

(horizontal) direction M = 2t = 12 between times t = and t = 5 + 1/2 (vertical direction). 

Down-pointing arrows denoting vacant sites have been omitted from the drawing. At time t = 

the even sublattice is filled and the odd sublattice empty. Particles are injected at site 1 after times 

t = and t = 4. At the boundary (site 12) particles get stuck at times t = 1 and t = 2 and are 

absorbed at times t = 0,3, 4, 5. 



With this convention Tl(cx, f3) acting on sites L and 1 corresponding to the vertex weights 
shown in fig. 3 is given by 

T L (a,(3) = 1 + a(s{ - a ± ) + P(st - r L ) + aP(sl - r L )(si - a ± ) 



( 



\ 



(A4) 



I- a (3(1 -a) 
a 1 af3 P 

{l-a){\-(3) 
x a(l-P) 1 - P ) 
The transfer matrix T = T(a, P) acts parallel first on all even-odd pairs of sites (2j, 2 j + 1 ) 
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including the boundary pair (L,l), then on all odd-even pairs. Thus in the first half time 
step T even shifts particles from the even sublattice to the odd sublattice (so far it was not 
occupied) and then, in the second half step, T odd moves particles from the odd sublattice to 
the even sublattice again. As a result, we expect an asymmetry in the average occupation 
of the even and odd sublattice which is related to the particle current. In a model with 
transfer matrix T = x odd T even the asymmetry will be reversed, but there will be no essential 
difference in the physical properties of these two systems. 

A possible configuration of particles in a 12x12 lattice is shown fig. 4. Note that the 
presence of particles at site x = 11 and times t — 2,3 imply the existence of particles on the 
left edge of their light cones as long as they move in a region where the even sublattice is 
fully occupied, i.e. they move with velocity of light (two lattice units per time step) to the 
left. A particle on an even lattice site at some (integer) time t always implies the existence 
of a particle on the right edge of its light cone up to the boundary. 

The model has a particle hole symmetry. We denote by \xi, x 2 , ■ ■ ■ , x^) = Sxi s x 2 ■ ■ ■ s x N I ) 
the iV-particle state with particles on sites xi,...,xn (| ) is the state with all spins up 
corresponding to no particle). The parity operator P reflects particles with respect to the 
center of the chain located between sites x = L/2 and x = L/2+1 and the charge conjugation 
operator C = n£=i a j interchanges particles and holes and therefore turns a iV-particle state 
into a state with L — N particles. One finds 



In the bulk the particle current is conserved and can be obtained from the commutators 
of r 2x and r 2x _i with T. These relations play a crucial role in the construction of the steady 
state and the computation of the time-dependent correlation function. Defining the current 
operators j|™ n and j^t-i by 



(CP) T(a, (5) (CP) — T(f3, a) . 



(A5) 



J2x — T 2x<?2x+1 



(l<x< L/2 - 1) 



(A6) 



■odd 
I2x-1 



(1 - o- 2:E _20"2x-i)(l - r 2x T 2x+1 ) (2 < x < L/2 - 1) 
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a straightforward calculation yields (x ^ L/2): 

[T,T 2x -l] = T (t 2x -1 - (1 - (T 2x - 2 CT 2x -i)t 2x T 2x+1 ) 

-odd 'even 

~~ J2x-1 J2x-2 

(A7) 

[T, t 2x ] = T {p 2x - 2 <J 2x -\{\ - t 2x t 2x+1 ) - a 2x ) 

-even 'odd 

~~ J2x J2x~\ ■ 

Current conservation implies that the expectation values of the current operators j|^ en 
and do not depend on x, (j 2x en ) = (j 2x -i) — const = j. 

Note that the cases a, (3 = 0, 1 are trivial. If a = no particles are injected and the 
steady state is | ). If a = 1 then in each time step a particle is injected and therefore the 
even sublattice fully occupied. Particles on the odd sublattice are randomly distributed 
with average density 1 — j3. As discussed above they move with velocity of light everywhere. 
Therefore the connected time-dependent two-point function on the odd sublattice is except 
on the left edge of the (forward) light cone. 
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